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Abstract 


The  boundary  value  problem  of  heat  conduction  In  a  three  dimensional, 
laminated  plate  is  approximated  by  a  hierarchy  of  two  dimensional  models. 
Computable  a-posterlori  indicators  and  estimators  of  the  modelling  error  in 
various  norms  are  derived  and  their  local  spectral  and  asymptotic  exactness  is 
proved.  Sharp  estimates  for  their  effectivlty  Indices  are  also  obtained. 


AMS  (MOS)  Subject  Classification: 

Primary:  65N15  Secondary:  65N30 


INTRODUCTION 


The  modelling  of  the  elastic  behavior  of  thin  objects  has  a  long  history. 
The  main  Idea  is  to  replace  the  particular  three  dimensional  problem  by  a  two 
or  one  dimensional  one,  which  is  easier  to  solve.  Such  approaches  were 
already  proposed  in  the  first  half  of  the  nineteenth  century  by  S.  Germain  [1] 
and  G.  Klrchhoff  [2].  Since  then  many  approaches  were  proposed.  For  some 
surveys  we  refer  for  example  to  [3]  and  [4].  The  derivation  of  these  models 
is  based  on  physical  considerations,  a  mathematical  analysis  of  various 
degrees  of  rigor  or  on  the  asymptotic  analysis  of  the  three  dimensional 
problem  as  the  thickness  of  thfe  structure  tends  to  zero.  We  refer  to  [5]  and 
references  there  for  this  approach.  In  general  all  the  methodologies  can  be 
understood  as  the  application  of  a  dimensional  reduction  approach. 

This  approach  leads  to  an  approximate  solution  of  the  original  higher 
dimensional  problem.  Hence  an  error  estimate  is  needed.  There  are  presently 
various  a  priori  error  estimates  (see  e.g.  [6],  [7])  or  estimates  of 
asymptotic  character  (see  e.g.  [5])  available. 

Nevertheless  in  today’s  computational  environments  we  need 

a)  an  accurate  and  computable  a  posteriori  estimate  of  the  difference 
(error)  between  the  exact  solution  of  the  original  three  dimensional  problem 
and  the  dimensionally  reduced  one,  the  modelling  error,  and 

b)  a  procedure  which  leads  to  the  construction  of  a  hierarchy  of 
dimensionally  reduced  models  which  allow  to  solve  the  original  three 
dimensional  problem  with  a  prescribed  given  tolerance  or  accuracy  and  this 
procedure  has  to  be  adaptive  (we  remark  that  in  contrast  to  the  classical 
approaches,  the  adaptive  approach  leads  to  models  which  are  not  uniform 
through  the  entire  domain). 

As  in  today’ 8  adaptive  finite  element  approaches,  the  fundamental  part  of 
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the  adaptive  procedure  are  a  posteriori  error  estimates  based  on  local  Indi¬ 
cators  which  should  be  of  high  quality.  By  this  we  mean  that  the  estimator 
and  the  indicator  have  to  be  robust,  i.e.  their  effectivity  index  should  be 
reasonably  well  bounded  from  below  and  above  for  a  large  class  of  solutions 
and  should  be  asymptotically  exact  for  more  restrictive  classes  of  solutions. 

The  present  paper  addresses  these  quest icns  for  the  heat  conduction 
problem  on  a  thin  domain  when  the  material  is  homogeneous  or  laminated.  It 
gives  a  computable  a-posteriori  estimate  for  the  modelling  error  measured  in 
the  (weighted)  energy  and  L^-norm.  The  indicators  are  local  and  hence  very 
well  suited  for  adaptive  approaches.  Upper  and  lower  bounds  of  their 
effectivity  indices  are  also  obtained.  The  adaptive  procedures  based  on  this 
approach  will  be  discussed  elsewhere  (see  e.g.  [8]).  Let  us  now  outline  the 
contents  of  this  paper. 

In  Section  1  we  introduce  the  formulation  of  the  problem  and  the  main 
notations.  In  Section  2  we  Introduce  the  hierarchic  models  and  some  of  their 
basic  properties.  Section  3  addresses  some  abstract  functional  analytical 
results  which  will  be  employed  later.  Section  4  introduces  the  a-posterlori 
estimator  and  proves  its  basic  properties,  especially  the  upper  and  lower 
bound  for  its  effectivity  index.  Section  5  analyses  the  asymptotic  exactness 
of  the  estimator  as  the  thickness  of  the  plate  d— >0  and  Section  6  analyses 
the  spectral  asymptotic  exactness  when  the  degree  of  the  model  increases. 
Sections  4,  5  and  6  address  the  estimator  for  the  modelling  error  measured  in 
a  weighted  energy  norm  where  the  weight  is  exponential.  Section  7  addresses 
the  error  estimate  for  the  measure  of  the  modelling  error.  Section  8 

generalizes  the  estimator  to  laminated  materials  and  the  final  Section  9 
presents  a  simple  numerical  example  to  illustrate  the  sharpness  of  our 
estimates. 
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1.  Notation  and  Problea  Formulation. 

2 

By  w  c  R  we  denote  a  bounded  domain  with  a  piecewise  smooth  boundary 
7  satisfying  the  cone  condition.  With  a  positive  thickness  parameter  d  and 
<i>  we  associate  the  three  dimensional  domain 


Q  ■  u  x  (-d/2,  d/2) 

with  lateral  boundary 

T  ■  7  x  (~d/2,  d/2) 

and  the  faces 

R±  =  {(Xj.Xg.Xg)  |  (Xj.Xg)  €  (<>  ,  Xg  =  ± d/2) . 

In  (2  we  consider  a  heat  conduction  problea  with  prescribed  heat  fluxes  f J 
on  the  faces,  i.e. 

Lu  =  0  in  Q, 

(1.1)  u  =  0  on  T, 


D  u  *  f' 
n 


on  R+. 


where  the  operator  L  is  (in  the  sense  of  distributions)  given  by 

■  (•£)  (»"•■)■ 


(1.2) 


-Lu  = 


aX3 


T 

where  7x  =  ,  g|-j  ,  x  *  (Xj.Xg);  a(»),  b(«)  c  L*(-l,l)  are  even 


functions  independent  of  d  and  satisfy 


(1.3)  0  <  A  S  a(z)  ,  0  <  B  i  b(z)  . 

The  matrix-function  C(x)  is  symmetric  and  uniformly  positive  definite,  i.e. 
there  exist  constants  0  <  C  S  C  <  »  so  that 

(1.4)  CISI2  S  CTC(x)C  S  CICI2 
2  —  €0 

for  all  £  e  R  ,  x  c  u;  C(x)  has  C  coefficients  and  the  boundary 
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operator  Du  is  the  (distributional)  exterior  conoraal  derivative  on  R.. 

n  ± 


To  cast  (1.1)  into  the  weak  fora  we  introduce  the  Sobolev  space 

(1.5)  H  :■  -^u  €  H* (Q)  | trace  u  *  0  on 

and  define  the  bilinear  fora  B(*,*)  :  HxH  — > F  and  the  functional  F(»)  : 

H — >R  by 

(1.6)  B(u.v)  -  %  %  *  ■>  £r)  <V,Tc(It>v}dItidx2dx3 

and 

(1.7)  F(v)  - 
respectively. 

Then  the  weak  fora  of  (1.1)  reads:  Find  u  €  H  such  that 

(1.8)  B(u,v)  »  F(v)  Vv  e  H. 

Under  the  assuaptions  (1.3),  (1.4)  there  exists  a  unique  weak  solution 
of  (1.8)  provided  that 

(1.9)  f+,  f"  c  L2(u) 

(this  assuaption  could  be  weakened,  but  is  sufficient  for  our  purpose). 

2.  Hierarchical  Modelling. 

We  will  approxlaate  the  boundary  value  problea  ( 1. S)-( 1. 9)  by  a  sequence 
of  two  diaenslonal  probleas  on  u,  the  hierarchy  of  plate  models,  which  we 
now  define. 

Denote  by 

(2.1)  f  -  (WjJwj  £  u  ,  1  £  i  £  n> 

a  collection  of  n  doaalns  with  piecewise  saooth  boundaries  du^  such  that 


t 

w 


f^Xj.Xg)  V  (Xj.Xg.l)  +  ^(Xj.Xg) 


V  (x 


l,x2 


’'l)}dxldx2’ 
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n  __ 

w.  n  <i>  *  f  if  i  *  J  and  u  -  U  v.  If  could  be,  for  example,  a 
J  i*l  1 

triangulation  of  u).  For  a  vector  q  of  nonnegative  integers 

(2.2)  q  =  (qj.-.-.q^  ,  qt  iO 
and  a  dense  sequence 

(2.3)  {Vz)}  C  h1(_1>1) 

J=0 

of  linearly  independent  functions  we  define 

q‘  Ov 

(2.4)  S(J>.,)  {u  4  H|u|Ui  -  £  V  #jP].  4 

J=0 

Then  S(T,  q)  c  H  and  the  (?,  q) -plate  model  is  the  boundary  value  problem: 
Find  u(y,q)  €  S(?,q)  such  that 

(2.5)  B(u(?,q),v)  *  F(v)  W  6  S(?,q), 


i.e.  u(?,q)  is  the  (energy)  projection  of  the  weak  solution  u  onto 
S(?,q).  Hence  (2.5)  constitutes  an  elliptic  boundary  value  problem  on  u 

(  4  ) 

for  the  coefficient  functions  Uj  (Xj.Xg)  in  (2.4). 

The  selection  of  the  functions  in  (2.3)  completely  determines  the 

(P, q)-model  and  has  been  investigated  in  19),  from  where  we  quote  the 


following  results.  Define  ^2j(z)  *  ^2J^~Z^’  ^  *  °*1'*-* 

(2.6)  a(z)^(z)v'(z)dz  -  0  , 

J-l 

(2.7)  a(z)^j(z)v'  (z)dz  +  b(z)*2j_2(z)v(z)dz  1 

■'-1  *-l 


recursively 


«,(*> 


for  all  v  c  H1(-l,l)  ,  J  e  M,  where 
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aj(v) 


(vd)  ♦  v(— 1 )  if  j  -  i 

0  else 


and  ^2J+l(z)  "  ■*2J+l(”z)’  J  ■  0,1,2,...  by 
(2.8)  a(z)^'2j+1(z)v/ (z)dz  +  b(z)#2j_j 

-1  J-t 


(z)v(z)dz 


for  all  v  «  H  (-1,1).  J  (M,  where 


5  (v)  -  {v(1>  -  v<-1) 
J  [  0 


If  J  -  1 
else 


(2.9) 


-  0. 


Remark  2  1 .  It  is  not  hard  to  see  that  (2.6)-(2.9)  determines  the  sequence 
uniquely  (the  nonuniqueness  In  the  solutions  of  (2.7).  (2.8)  is  taken  care  of 
by  requiring  the  compatibility  condition  in  the  subsequent  step).  Moreover, 
it  was  shown  in  [9]  that 


(2. 10) 


{*j(z)>J»0  lS  dense  ln  H (-1,1). 


[.  If  a(z)  and  b(z)  are  constant.  #j(z)  Is  a  polynomial. 


Table  1  lists  the  first  four 


Vz} 


(3z*  -  1 )/6 
(z3  -  3z)/6 
(15z4  -  30z2  ♦  7)/360 


Table  1.  The  first  ^(z)  for  a  ■  b  »  1,  0  S  J  S  4. 
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If  a(z),  L(z)  are  piecewise  constant,  ^j(z)  is  *  piecewise  polynomial. 
This  the  situation  for  sandwich  materials. 

With  the  choice  (2.6)-(2.9)  of  ^j(z),  the  modelling  error 

(2.11)  e(P,q)  :*  u  -  u(P,q) 

Is,  for  f  ■  {«>  and  q  ■  N,  of  optimal  asymptotic  order  as  d— >0, 
provided  the  data  f+,  f  are  sufficiently  regular  In  u  and  satisfy 
certain  compatibility  conditions  on  the  edges  dr  *  y  x  {±d>  which  ensure 
the  absence  of  boundary  layers  (see  [10],  for  example).  We  emphasize  at  this 
point  that  due  to  (2.10)  the  error  e(!P,q)  will  also  tend  to  zero  for  fixed 
d  >  0,  if  min{q^>  — »»  In  contrast  to  the  error  in  models  obtained  by 
asymptotic  analysis. 

In  the  following  sections  we  will  derive  computable  a  posteriori  error 
estimators  for  the  modelling  error  measured  in  the  energy  norm 

(2.12)  |e(?,q)BE(Q)  “  <B(e,e) )l/2 

In  terms  of  the  residual  data  on  the  faces  R±.  The  following  property  of 
e(P,q)  will  prove  to  be  Important. 


(2.13) 


L.  For  every  (?,  q)  we  have 


I  d/2 

b[?]e<XfVX3)dX3 

/o 


a.e.  (Xj.Xg)  c  w. 


Proof.  It  follows  from  (2.6)  that  0q(z)  “  const.  Due  to  (2.5)  we  have 

B(e(?,q),  v)  »  0  Vv  e  S(P,q) 

and,  since  X(x)^Q|-~j  «  S(P,q)  for  all  (P,q),  we  find  with  Fublnl’s 
theorem 
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r 

(2.14) 

0  "  VxX(Xi*x2)*C(x1-^)  7X 

bme(xi>x2'x3)dx3dxidx2 

J« 

-d/2 

for  all 

X  €  ft^u)  *  (u  €  H*(«)  |  trace 

u  on  r  =  0>. 

Let 

* 

*  bl- 

*3l 

d-Je(Xl*X2>X3)dX3  * 

J-d/2 

Since  e 

€  H,  *  €  ft1 (Q) .  Using  X  =  * 

in  (2. 14)  we  conclude  that  ^  *  0 

which  was  to  be  proven.  □ 

We  shall  derive  in  Section  4  computable  guaranteed  upper  estimates  for 
the  modelling  error  (2.12)  which  are  asymptotically  exact.  Moreover,  our 
estimators  also  give  information  about  the  local  contributions  from  to 

^eH£(Q)  As  a  tool  we  shall  use  certain  exponentially  weighted  spaces,  which 
we  analyze  next. 

Remark  2.3.  Here  and  in  what  follows  we  assume  that  the  elliptic  system  (2.5) 
of  differential  equations  for  the  unknown  coefficient  functions  Uj1^  is 
solved  exactly.  However,  usually  only  an  approximate  solution  can  be 
obtained.  Our  results  remain  nevertheless  valid,  if  the  approximate  solution 
has  a  sufficiently  high  accuracy.  In  computational  practice  one  usually  works 
with  a  finite  element  approximation  of  In  this  case  the  a  posteriori 

estimation  of  modelling  error  can  be  used  to  determine  the  desirable  accuracy 
of  the  finite  element  approximations  of  Uj^.  □ 
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3.  Some  abstract  results. 

Let  Hj,  Hg  be  two  reflexive  Banach  spaces  furnished  with  the  noras 
1  •  H j  and  I’llg*  respectively.  Further,  let  B(u,v)  be  a  bilinear  fora 
defined  on  Hj  x  Hg.  We  will  call  the  bilinear  fora  (C, ^-regular  if  there 
exist  constants  0  <  C,  jr  <  «  so  that 

(3.1)  |B(u,v)|  S  C|u|j  |v|2. 

(3.2)  inf  sup  |B(u,v)|  t  r, 

lulj-l  flv|)2=l 

(3.3)  for  any  v  *  0,  v  €  H_,  sup  |B(u,v)|  >  0  . 

Rulj-1 

Bilinear  foras  satisfying  (3.1)  -  (3.3)  have  the  following  properties. 

a)  Let  f  €  OL,)'  (1*e*  Is  a  bounded,  linear  functional  on  Hg), 
then  there  exists  exactly  one  u  c  such  that 

B(u,v)  ■  f(v),  Vv  e  1^. 

b)  If 

(3.4)  sup  |B(u,v)|  S  A, 

Iv#2-1 

then 

|u| ,  s  -  . 

1  r 

Let  us  consider  now  soae  special  cases  which  will  be  laportant  later. 

Let  0  <  pfoj'Xg)  «  W1,w(«)  denote  a  (strictly  positive)  weight 
function  in  w.  We  define 

H  ■  (e  c  H  I  e  satisfies  (2.13)} 

9 

and  furnish  with  the  weighted  energy  norm  defined  by 

(3.6)  |.|J  -  f  f2<vV  {afr]fe)  * b  (-a2)  VTcv}dxidx2d*3  • 

JQ 

The  following  Leaaa  will  be  used  repeatedly. 
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Leaaa  3. 1.  Assume  that  u  c  H^.  Then 

(3.7)  |  ^b^j^dx^dxj  s  A2d2  |  *2afr](l^)  *'1'*V*3 

<rx( -d/2,  d/2)  <rx  ( -d/2 ,  d/2 ) 

for  all  open  subsets  <r  £  ».  Here  A  is  given  by 


inf 

*€^(-1.1) 


and  the  lnflaua  is  taken  over  all 


1>  €  H1  (— 1. 1 )  a  |^|J  b(z)^(z)dz 


Proof. 

Assuae  that  u  e  C*(G)  n 

H  .  Then  we  have  for  all 
P 

x  e  u  the  bound 

rd/2  ^x_^  , 

> 

(3.8) 

f  b(^)lul  dXj 

% 

*3 

-d/2 

J-d/2  ^ 

by  the  definition  of  A  and  a  scaling  arguaent.  Multiplying  both  sides  of 
2 

(3.8)  by  9  and  integrating  over  <r  we  get  (3.7)  for  u  and  a  density 
arguaent  coapletes  the  proof.  □ 

2 

Remark  3. 1.  For  a  *  b  *  1  we  find  A  *  -  .  If  ^  is  syametrlc  in  z 
then  A  ■  — . 

it 


1.®/ 


Theorea  3.1.  Let  0  <  p(Xj,x2)  e  W  ’  («)  and  assuae  that 


Q  I  shXr~l- 


l»i,2 


<  oo. 


L  («) 


Define  Hj  ■  H^,  Hg  ■  Hj^.  Then  the  bi linear  fora  (1.6)  is  (1,7)  regular 
on  Hj  x  Hg  with 
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(3.9)  yir0  £l  -  d AQ  / C/2  j  £l  +  d  /ST  AQ^l  +  d*4cA<lj"J 


Proof: 


1)  |B(u,v)|  5  |uS^  lv| follows  immediately  from  Schwarz’  Inequality. 

2 

2)  Let  us  show  (3.2).  For  u  c  define  v  «  p  u.  Then  v  e 


_  2_  _  ,  2*  flv  2  5u 

V  ■  *  V*  *  UV*  >•  irg  -  * 

Hence,  denoting  the  volume  element  by  dx,  we  find 


(3.10) 


B(u,u)  -  |u|~  + 


{b(-^]uVx(p2)TC(x)VxuJ.  dx, 

»  O 


and  we  estimate  for  every  c  >  0 

|  J  {b(ir)“  y*2>Tc<*> v} 


J  Qv  C/2 


■  F  f  b[^ju2f»2<Jx  +  c  f  P2b |^jp) 7xuTC ( x ) Vxu  dx  •. 


Utilizing  Lemma  3.1  with  <r  *  w,  we  arrive  at 


|  {b(-?]uV*2,Tc(x)V*} 


4  ,2b^r]vTc,J',Vd1'} 


Selecting  c  ■  Ad  yields 


B(u, v)  2  (1  -  d  AQ 


/ C/2)  |u|2 


Further, 


Ivfl 


2 

l/f> 


|u|2  +  f  f“2b[^j/2f2u7x(f2)TC(x)7xu  +  u2Vx(f»2)TC(x)Vx(f2)|  dx 

JQ  J 

S  |u|2  +  dKQ/F  |u|2  +  2d2A2C  f  f2a(^j^g-]2dx 
£  +  d AQ^l  +  d AqJJ  |u|2 


from  where  get  (3.9). 


□ 


Remark  3.2.  We  observe  that 


Q  *  2  max 

1*1,2 


3x^/  9 


L°*(«) 


and,  from  (3.9),  we  can  select  p  in  particular  so  that  a 
and  get 


d  <  1 


(3.11) 


r  *  (1  -  a)  (1  +  2a(l  ♦  2a))“J/Z 


4.  A  posteriori  estimation  of  the  modelling  error. 

In  this  section  we  assume  that  the  (f.q)-model  (2.5)  has  uniform  order 
q  and  that  Its  exact  solution  u(?,q)  Is  known.  We  will  be  Interested  in 
computable  estimators  for  (2.12),  the  modelling  error  in  energy  norm.  To 
avoid  obscuring  the  main  ideas  behind  technicalities,  we  assume  throughout 
this  section 

(4. 1)  a(z)  -  b(z)  -  1  ,  C(x)  -  £  J  °  J  , 

l.e.  L  in  (1.1)  is  the  Laplace  operator  and  *  8/8 n  is  the  outside 
normal  derivative.  All  results  apply  with  minor  modifications  In  the  proofs 


12 


which  we  will  present  in  Section  8  also  to  (1.1).  It  is  convenient  to  write 
u  =  ul  +  where 

(4.2)  UjfXj.Xg.Xg)  «  -UjUj.Xg,  -  Xg),  UgtXj.Xg.Xg)  «  UgUj.Xg,  ~  Xg  )  . 

The  Uj  satisfy  :  c  Hj  such  that 

(4.3)  B(uifv)  *  F^v)  Vv  e  Hr  i  =  1.2 
where 

Fj(v)  *  f1(x1,x2)(v(x1,x2,d/2)  -  v(x1,x2.  -d/2))  dXjdXg  , 

F2(v)  ■  f2(x1>x2)(v(x1,x2,d/2)  +  vtXj.Xg,  -d/2))  dXjdXg  , 

Jw 

(4.4)  fjUj.Xg)  -  i  (f+  -  f')(x1,x2).f2(x1.x2)  -  \  (f+  ♦  f^JtXj.Xg)  . 

Hj  »  (u  €  H  |  u  is  antisymmetric  (symmetric)  in  Xg  for  i  *  l(i«2)>. 
Obviously,  the  spaces  Hj  and  Hg  are  orthogonal  in  energy,  i.e. 

(4.4)  B(u, v)  »  0  V  u  6  Hj,  Vv  6  ftj  , 

and  u(?,q)  ■  u^f.q)  +  UgfP.q),  each  of  which  can  be  obtained  by  energy 
projection  of  u^  onto 

(4.5)  Sj(f,q)  :«  S(?,q)  nHj,  i  ■  1,2. 

Further,  from  (4.4)  we  get  also 

(4.6)  |e(f,q)|g(Qj  *  lei^*<l)l£(Q)  +  le2^‘<I^E(Q)  * 

where  e^f.q)  *  -  Uj(f,q),  1  *  1,2. 

Since  S(f*,q)  in  (2.4)  depends  only  on  spanl^jU)}^  ■  n^t-l.l)  (see 
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Remark  2.2),  we  will  assume  below  for  convenience  that 


(4.7) 


#j(z)  *  Lj(z)  , 


where  Lj  denotes  the  Jth  Legendre  polynomial  on  (-1,1). 


All  our  a-posteriorl  estimators  8  for  (4.6)  are  of  the  form 

(4.8)  g(u. 


^(P.q))  *  iDjfXj.Xg)  ^dXjCb^ 


,  1  *  1,2. 


Here  ^(Xj.Xg)  Is  called  Indicator  function. 

Let  |o|  be  any  norm  on  H  and  8  In  (4.8)  an  a-posterlorl  error 
estimator  for  ge(9>, q)|.  Then  we  define  the  effectlvlty  index  8 
corresponding  to  8  and  |og  by 


(4.9) 


«  S(u(f,q)) 

®  •  leTFTqTir 


Vfe  say  that  8  is  a  guaranteed  upper  estimator.  If  6M  for  all  u.  The 
estimator  8  Is  (k ^K^-proper  with  respect  to  a  class  T  of  data.  If 

0  <  ic1  S  8  S  k2  <  «  Vf  «  T. 

Further,  8  is  asymptotically  exact  on  T  If 

(4.10) '  8  - »  1  as  d  - »  0+  Vf  €  T, 

and  8  is  spectrally  exact  on  T  if 

(4. 10) *  8  - »  1  as  q  - >  »  Vf  e  T. 

Finally,  8  is  locally  asymptotically  (resp.  spectrally)  exact  on  T,  if 

(4.10) '  (resp.  (4.10)*)  hold  with  the  norm  |«|  defined  in  (3.6)  where  the 

9 

weight  function  ^(Xj.Xg)  is  given  by 

^(Xj.Xg)  :■  exp  “p^lxi-x2l  ♦  IXg-x^lj,  0  <  p  S  1,  0  S  cc  <  l/(Av£E~' ) 
and  (x*,x^)  c  u  is  arbitrary. 
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We  begin  the  analysis  of  the  estimator  S  for  the  case  9  •  {«>  and  the 

2 

energy  norm  and  consider  the  L  -norm  later  in  Section  7.  Whenever  the  order 
of  the  model  is  uniform  throughout  «,  we  shall  omit  the  index  9.  Due  to 
(4.6)  we  can  derive  the  indicator  functions  ijj,  i  =  1,2,  separately.  We 
start  by  observing  that,  due  to  Theorem  2.1,  the  errors  e^Cq)  e  H^,  i  *  1,2, 
defined  in  (4.6)  satisfy 

(4.12)  B(ei(q) , v)  *  R^v)  Vv  e  H 
and 

(4.13)  B(eltv)  *  0  Vv  €  S^q) 

where 

Rj(v)  *  ri(x1,x2)  (vfXj.Xg.d/B)  ±  v(x1,x2>-d/2))dx1dx2 

+  v(x1,x2,x3)  AUj (q)dx1dx2dx3,  i  «  1,2 

•'ll 

and  -,+  correspond  to  i  «  1,2,  respectively.  Here 

Su. (q) 

(4.14)  ri(x1,x2)  =  fi(x1,x2) - ^ —  (Xj.Xg.cL^),  i  *  1,2. 

Remark  4. 1.  For  fj  e  slO,  we  have 

(4.  IS)  r,  *  H,lln(s-I*C,(U) 

where  c  >  0  is  determined  by  the  maximal  regularity  of  the  solution  to  the 
Dirichlet  problem  for  in  u  (e.g.  if  u  is  a  slit  domain,  0  S  c  <  1/2). 

The  unknown  coefficient  functions  in  (2.4)  satisfy  the  elliptic  system 

-|AAxU  +  |BU*  cV  -  cY  in  u 

U  »  0  on  8u. 

The  matrices  A  and  B  are  Independent  of  d  and  given  by 
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and 


•u  ■  f  * 
J-1 


v*’  B1J 


I  *[*y iz 


c*  ■  <^0(±1),...,  *q(±l)>T 


We  calculate  next  a  slaplified  expression  for  R^(v)  which  we  will  use 
repeatedly  below. 


Lemma  4. 1.  Let  i  “  1,  q  =  2m+l  or  i  ■  2,  q  »  2m  then 

R. 


(4.16) 


j(v)  *  TjfXj.Xg)  |v(x1,x2>d/2)  ±  vtXj.Xg.-d^) 

-  [  V(X1,X2.X3)  g|- 

-d/2 


Proof:  Let  i  *  1,  q  «  2a+l.  Then 

2m+l 


(4.17) 


iUj(q)  -  £  LJ  £r] 

J=0 


,“1, 


for  some  Ajj  €  H  («).  To  determine  Ajj,  we  use 

Rj(v)  *  0  Vv  €  Sj(q)  u  Hg. 

We  select  v  *  V(Xj,Xg)  L^  *  Hg  wlth  arbitrary  V  e  ft1(«)  and  get 

Ajj  «  0  for  even  J.  For  J  *  2k+l,  0  S  k  S  a,  we  find 

0  -  f  VIXj.Xj,)^,  *  A1Zkti  |  dXjdXj. 

-1 

Since  we  get  V  e  ft^w)  is  arbitrary,  we  get 


A1.2k.l-  *  I  (2(2k*1)*1>  rl 
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(note  that  due  to  Remark  4.1,  Aj  2k+1  e  L  («)). 


Hence 


*!,(,)  -  -  |  £  C4k+3)  1^,,^) 

k-0 

'  -  rl(Vx2>  [4..2 (-?))• 


For  i  *  2  and  q  -  2m,  one  proceeds  analogously. 


We  derive  next  the  estimator  6.  We  start  with  the  observation  that  from 
Theorems  2.1,  3.1  we  have  the  bound 

IBUj-Ujtqhv)  | 


(4.18) 


7n  lu.  -  u  (q)|  S  sup 


|v| 


1/p 


where  rQ  Is  as  in  (3.9).  With  (4.12)  thus 


-  ,  (R.  (v))‘ 

y^»e.(q)|j  S  sup  - * — 


*  °*'v»i /,  |v|T/» 


S  sup 
0*«v| 


(Rt(v))‘ 


U*  I/’2 (%]  (W* 3,dx 

S  sup  f  f  pr.  *, [v]  dx.dx-l 

0*lv|.  LJW  2J 


1  »  1.2 


where 


(*iIvl]{xrx2) 


(4.19) 


vlXj.Xj.d^)  ?  Wx.^.-d/Z)  -  J^vCXj.Xj.Xj)  ^ 


(C(%h 


1/2 


Hence  we  obtain,  using  Jensen’s  inequality,  that 
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(4.20) 


s  }  ^  (*1tvl)2dXl<lx2 


where  the  supremum  is  taken  over 


(4.21)  M  :«  L 


2(«,H1 (-d/2, d/2) )  n  |v|  v  dXg  -  0  a.e.  (Xj.Xg)  €  »| 

—A  /O 


(see  (111  for  the  definition  of  anisotropic  Sobolev  spaces). 

Since  is  strictly  concave  and  upper  seal continuous  on  M,  there 

exists  a  (unique)  aaximizing  eleaent  £  M  which  satisfies  the 
Euler-Lagrange  equations 


_  _d_ 
dx3 

(viRl 

)  ln 

(-  !•!)• 

I 

.(♦> 

if 

i*i. 

l±d/2 

1  ±1 

if 

i  «  2. 

Hence  we  find  that  Vj  is  Independent  of  (Xj.Xg)  and  given  by 


(4.22) 

and 

(4.23) 


,2X3*1  ^2X3* 

*  *  V2IT]  ~  LqL"d~]  q  *  o 

i  2  2q+3  *  q  ’ 


(4.23)  (*ilvi1) 

Referring  to  (4.20),  we  have  proved 


•  l2  d 
1  «  :  dC, 

i  I  2q+3  iq 


Theorem  4. 1.  Assume  that  fj  in  (4.3)  is  square  integrable  over  «.  Then 
the  error  (q) 8e^qj  for  the  hierarchical  model  of  uniform  order  q  (l.e. 
f  «  {«>  and  odd  q  1  1  for  every  i  *  1,  q  even  fo~  i  *  2)  can  be 
estimated  by 


(4.24) 


s  shs  dxidx2 

«  ik 


where  f  is  as  in  Theorem  3. 1  and 
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ri(xrx2)  "  fi(xrx2)  " 


(Xj.Xg.Xg),  i  «  1,2. 


dx3 

Based  on  (4.24)  we  define  the  indicator  functions 

(4.25)  Hi(j  (Xj.Xg)  »  J 2^3  ^(Xj.Xg)  r^lXyX^),  i  *  1,2 

and  the  estimates  6(u^(q))  defined  in  (4.8)  are,  according  to  (4.24), 
guaranteed  upper  estiaates  for  (e^q)!^^,  since  p  ■  1  iaplies  tq  *  1  ln 
(3.9). 

Remark  4.2.  We  emphasize  that  is  very  easy  to  compute,  especially  for 

low  order  models.  We  find  in  particular  for  i  *  2,  q  *  0  that 

1*20*xrx2*  *  3  WXj.XgJfgtXj.Xg)) 

Selecting  p  ■  1  Implies  Q>0  in  Theorem  3.1,  whence  we  obtain  yQ  *  1  in 
(3.9).  Thus  (4.24)  yields  with  r2  *  f2  *  (f+  ♦  f  )/2  the  estimate 

(4.26)  '«2<0'l|(Q)  S  12  ,fV''^(u)- 

In  the  subsequent  sections  we  will  demonstrate  that  the  estimators  £(Uj(q)) 
based  on  (4.25)  are  asymptotically  and  spectrally  exact. 


5.  Asymptotic  exactness  of  the  error  estimator. 


Before  demonstrating  the  asymptotic  exactness  of  S,  we  introduce  some 
notation.  Throughout,  p  will  denote  the  exponential  weight  function  (4.11) 
Further 

(5.1)  l<,  -  {  I'frlV  . 

Finally,  we  Introduce  the  class  of  data 

(5.2)  T  :«  {fleither  r.(f)  *  0  or  Ir.l,  /|r, |n  £  0  <  •}. 

p  i  i  i»P  l  u»P 

The  main  result  on  asymptotic  exactness  is 


Theorem  5. 1.  Let  0^,  1  »  1,2  denote  the  effectivity  indices  (4.9)  with 
respect  to  the  weighted  energy  norm  (3.6).  Assume  further  that  f  »  {«},  l.e 
the  model  order  is  uniform.  Then  for  i  *  1,2  holds: 

1°  If  fj  e  L2(u)  we  have  with  Aj  »  |  ,  Ag  «  ^  ,  that 

(5.3)  Sj  i.  *n  :«  Jl-Aj  qJ  ^1  ♦  dVS^AjQd  +  d  V^A^I  ,  1  «  1,2  . 


2.  If  fj  c  Tp  then 


(5.4)  0J  S  x12  [l  ♦  |  d2Dq  (p2  +  Q2)J*  ,  i  «  1,2 


where  is  given  by 


(5.5) 


q  (2q+3)Z  -  4 


Moreover,  if  pal,  the  factor  3/2  In  can  be  replaced  by  1/2  and 

Q  ■  0. 
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Proof:  1.  The  bound  (5.3)  follows  immediately  fro*  Theorems  4.1  and  3.1,  if 


we  note  that  Aj  *  2/k  and  *  1/w  in  Lena  3. 1  with  a  *  b  *  1  since 

the  infimum  there  is  taken  only  over  odd,  resp.  even  #  €  H1 (—1, 1 ) . 

2.°  To  show  (5.4),  we  select  in  (4.12) 

—  2  •  2 
v  *  vp  *  v^Xg)  ritxl,x2)f> 

with  v^  as  in  (4.22)  and  get  with  (4.23)  that 

(5.6)  R^v)  -  dCq  |  r2p2  dx^  -  Bte^ql.v)  5  ge^  gvg^  . 


Since 


|Vv|2  S  p2|3f2|yxv|2  +  6 | v |  IVxPl2} 


we  find 


f  »“2{l  V'2  *  (%]  }  <b'ldx2dx3 


j  p17, 


V I2  ♦  6|v|2  | 


V2  *  »2  (%)  }  dxld»2dx3 


Since 


f  iv|i2[3»2  i  Vi  I2  *  6IV|Z  lri|2]  *  (et)  *2|^|Z } ‘‘VV’V 

JQ  J 

r  -  *. 


jd/2 

|  (v;)2dX3-|d3CqDq  , 

j  -<(  /O 


where  D  is 
<1 


in  (5.5),  we  get  with  |Vxf|2  *  Q2f2/2 
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The  local  asyaptotlc  exactness  not  only  ensures  that  the  Indicator 
functions  in  (4.25)  give  a  good  estimate  of  the  global  modelling  error 

A# 

in  energy  norm,  but  also  for  subdomains  u  of  u  that 
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*iq(p;w) 


l 

5 


r^Xj.Xg)  dx^ 

3^;«»)<dp 

Is  an  asymptotically  exact  measure  for  the  local  contribution  to  the 
modelling  error  at  u. 

It  further  guarantees  that  a  local  increase  of  the  model  order  in  u 
will  reduce  the  error  in  u  while  leaving  the  error  elsewhere  in  u\u 
unchanged  —  a  feature  typically  not  found  for  elliptic  equations  and  a 
consequence  of  the  fact  that  Q  is  a  thin  domain.  These  observations  are  the 
basis  for  the  adaptive  selection  of  the  model  orders  on  subdomains  of  u  [8]. 

6.  Spectral  exactness  of  the  error  estimator. 

Our  purpose  in  the  present  section  is  to  show  the  spectral  exactness  of 
the  errcr  estimator  S,  i.e.  that  k^,  in  Theorem  5.1  tend  to  one  as 

q— »»  at  fixed  d  >  0.  While  this  is  not  hard  to  establish  for  in 

(5.4),  the  corresponding  proof  for  »c^  requires  a  more  careful  analysis  of 
the  constant  A  in  Lemma  3.1. 

2 

We  denote  by  the  eigenfunctions  (orthonormal i zed  in  L  («)) 

of  the  eigenvalue  problem 

(6.1)  -Af>k  *  Akf>k  in  w  ,  p k  «  0  on  du  , 

k  *  1,2,3...  and  the  eigenvalues  are  enumerated  with  respect  to  increasing 
magnitude  and  counting  multiplicity.  We  collect  some  of  their  properties  that 
will  be  needed  later. 

Lemma  6. 1.  Assume  that  u  c  Rn,  n  ■  1,2  and  that  du  is  smooth  if  n  ■  2. 
Then  lT  >  0  is  a  simple  eigenvalue  and  the  corresponding  eigenfunction 


dist(x. 
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C  U 


fjtXj.Xg)  >0  in  «. 


2?  For  all  k  c  M,  and  all  (x^.x^) 


(6.2) 


■  n  5  C(«)0.  ,  0. 

If>l'c0(5)  k  k 


JX£/2  If  c*  c  R1 
[±+*  If  «  c  R2  " 


Proof.  Assertion  1*  is  well  known.  Let  us  prove  2°.  In  the  case  n  »  1 
the  assertion  follows  froa  the  explicitly  known  eigenfunctions.  We  consider 
therefore  n  *  2  and  claim  that 

Sf.l 

-  srl,  4  C1<“> >  0  • 

To  prove  it,  we  note  that  >  0  in  u  and  hence  the  function  u  * 
satisfies 


Au  *  Ajfj  2  0  in  u  ,  u  <  0  in  u,  u  *  0  on  du. 

Now  we  can  apply  the  aaxlaua  principle  (12,  Theorem  2.7]  and  find  that 

(xQ)  >  0  for  all  xQ  €  du  from  where  the  claim  follows. 

Now,  since  y  is  smooth,  there  exists  c_(«)  >  0  such  that  u  :* 

u  c 

{x  e  u|dlst(x,y)  >  c>  is  smooth,  too,  for  0  <  c  <  cQ.  Hence  we  can 


estimate,  possibly  after  reducing  cQ, 


and  get 


\>\ 


f>j(x)  2  i  Cj(u)e  on  wc 


LI..,  S  ^ 


1.  e. 
(6.3) 


Vu>  'VVL2 


(u) 


$ 


s  c3(«)  l»klL2(u)  -  c3 
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In  «r  :*  u\u  we  introduce  the  boundary-f i  1 1  ed  coordinates  (p, s) .  Then  we 
c  c 

have  fro a  Taylor's  formula  that 


|f>k(p,s)| 


IPjtp.s) |  8f 


2 

I Pa~  (0,s)  +  0(p  )| 

°P - -  S  C4(»)  sup 

2m  xer 


I (0,s)  +  0(p  )  | 


* 

ST  (x) 

9*1 

- 

k 

a r  (x) 

4 

,  c4(«)  |s»k  I  c4  |a,k  I 
s  (TT5T  I  an  I  s  c  cs<">  laiH 
1  1  *L  (r)  1  1 


S  C6(U)  ,’’klH2*2S(u) 


where  Cg  is  independent  of  k,  but  depends  on  5  >  0.  Hence 


(6.4) 


^  s  |a-\il2()  -  v-. 

co  (V 


Combining  (6.3)  and  (6.4)  completes  the  proof. 


□ 


With  the  eigenfunctions 

*i.  [r]  -  aln(Mi.  ¥\  ■  •‘i.  "  t-1  * 

[r] '  cos(^.  r]  •  "2. '  “ 

m  “  1,2,3,...  ,  the  eigenfunctions  for  -A  on  Q  with  homogeneous  boundary 
conditions  on  T  and  R±  which  satisfy  (2. 13)  are  given  by 

(6-7)  *ikm  WV^i.  (~d~ ]  ' 

with  the  corresponding  eigenvalues 

(6.8)  *lta  -  V44'2**?.'  1  *  >-2  k'  «• 


(6.5) 

and 

(6.6) 
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The  sequence  of  eigenfunctions  is  dense  in  the  space  H  defined  in 

(3.5).  Therefore  we  have  in  particular  for  the  modelling  errors  e^  the 
expansion 

(6.9)  «,  '  £  E,k»  #lta  •  1  * 

k.meN 

where  the  coefficients  can  be  detersined  froa  Lemma  4.1. 

Lemma  6.2. 

Eikm  =  (2  \  *  d  Mim)  Pim  plk 

with  and  as  In  (6.10)  and  (6.11)  below,  respectively. 

Proof.  We  recall  (4.12) 

B(erv)  -  Rj(v)  Vv  €  H  i  -  1.2, 

where,  with  (4.16)  and  integration  by  parts, 

d/2 

Vv)  *  f  ri(xrx2)  [  Lq+1  £r]  d*3d*ld,l2  • 

Ju  J-d/2  J 

Since 

B(W  ♦jin]  ‘  (I  *k  *  I  •‘?.]*lj*kl*» 

we  obtain  with 

(6.10)  Plk  «  I  ^(Xj.Xg)  ^(Xj.Xg)  dXjdXg 
and 

(6.11)  -  #'.(2)  Lq+1(z)  dz 

^-1 

the  assertion.  □ 
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Our  purpose  is  to  estimate  the  dependence  of  A  in  Lena  3. 1  on  the 
aodel  order  q  since 


A/2  rde 


(6. 12) 


inf  inf 


CJ?K  , 


(x1,x2)e«  e,  (e  }2  (AjfqJd)' 

-d/2  1  ^ 


where  the  inf i bub  is  taken  over  all  e^  of  the  fora  (6.9).  We  compute 
A/2 


L  J«i>2  -3 


’-d/2 


<•»>  i’ i  **.(*•#  »"(■•?  '±\ 

k,  l.n  ^ia  Mia 


,2  X,  ,-l  CB?.)2 


"u  ‘V 


iVi 


and 


d/2  £^2 

J-d/2  lax3J 


**3 


«■*«  ta  i 

k,  l,  a 


2  X,  ^-1  Ojj2 


la 


la 


"*i.) 


2  fk9l 


To  obtain  a  lower  bound  for  (6.12),  we  estimate  (6.13)  froa  above  and  (6.14) 

2  2 

from  below  (pointwise).  For  (6.14)  we  have  with  a  ■  1  +  d  A1/(lp1J) 

A/2  faen2 

f-d/2  W  ^  * 

I  {ip,/!*/*-2  -  2,h'  •‘’ii1  Z  (>  *  r  Vf 'V 


kk2 


r~>  /•  J  1  •j'l/'  rf2  -i ~1  >  (P?  ) 

Z  Kt]  Kj)  '»u«  <»«'  '*«■}  ^ 

c^i2  Mil  *lm  lMim' 


q  |2 
2 


d  -2 
“2  a 


fl  -  2  Z  f 

1  +  — 5^ 

1  fe  1 
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(6.15) 


q  ^2 
2 


f (ft 

=  |  a"2|p  I2  If  |2  (1  -  2f  -  f2)  V  — ^ 

.  “V 


where  we  assumed  that  p^  *  0  and  defined 


(6.16) 


*  :*  C(w) 


k£2  4^il 


d^-i  IPlkl 


0. 


IPill 


with  C(<<>)  and  as  In  (6.2). 

Consider  next  the  upper  bound  for  (6.13).  We  estimate  analogously  as 
before  and  get 

fd/2  _  ,g(j  »2 

(6.  .7)  [  (e,)2  <1x3  S  [|]  I*,I2IPUI2  a'2  (1  *  *)2 

J-d/2  m  lMini 


1*1,  q*2m+l  or  i  =  2,  q  =  2m. 


Now 


0 


im 


1.  Vi(z) 


r 


Vll5)d5dZ  *  *i.  Lq+1  d? 


Vi(C,d?  dZ 


z=±l 


where  we  used  that  (±1)  *  0  and  *  p2^  Hence  Oj^)2/ 2 


la 


are  the  Fourier  coefficients  of  the  antiderivative  of  L  i.e.  of 

q+1 
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Therefore 


(2q  +  3)_1(Lq+2(z)  -  Lq(z))  ‘ 


(6.18) 


y 


K(2q+3)"2  |L  £  -  L  »% 

q  2  qL2(-l,l) 


K 


and  analogously 


(6.19) 


( (2q+3;  -  4) (2q+3) 


(P?,)2  2 

— *  K|L  r  *  K 

,2  q+1  l2(.,  „  2q+3 


where  K  Is  a  constant  depending  on  the  normalization  of  the  ifi,  in  (6.5), 

IB 

(6.6)  (its  numerical  value  Is  immaterial  in  what  follows).  We  can  now  prove 
Theorem  6. 1. 

l.°  If  f  ,€  T  _  defined  In  (5.2)  with  0  *  0(q/d)1-c,  c  >  0,  and  0 
i  P 

independent  of  d  and  q,  then 

— » 1+  both  as  d  — *0+  ,  q— . 

2?  If,  moreover,  f ^  is  such  that 

«)  Pn  *  0 

0)  ♦(Tj,  u,d,q)  S  1^2^  -  1)(1  -  Dq  c)  for  some  e  >  0,  then 

Kjj— >1  both  as  d— »0+  ,  q— >«». 

Proof:  Assertion  1°  follows  immediately  with  the  definition  (5.4)  of 
and  the  assumption  on  0.  To  show  2°,  we  note  that  ic^  is  as  in  (5.3), 
however  now  with  Aj  determined  from  (6.12)  instead  of  Lemma  3.1.  Using 
(6.15),  (6.17)  in  (6.12)  we  find 
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1  4_  1-2 *-r  _m _ 

(Al(q)d)2  d2  iW  £  ,^,2,^^ 


and,  froa  (6.18),  (6.19)  with  *  -  (VS^-1)  -  5,  0  <  8  <  V§^  - 


1  ^  4_  6(2V2~,-a)  (2q+3)  -4 

(A^(q)d)2  d2  (V§^-8)2  2 

^  (1  +  v^)  . 

d  u 

q 


Hence 


Aj(q)  5  (1  +  V^")_1D  /« 


and,  using  6  «  (v^  -  1)  D1  c,  we  find  In  (5.3) 

q 


Aj(q)  1DJ*  ( (2q+3)2  -  4)' 


This  coapletes  the  proof. 


Therefore,  under  the  assumptions  Bade,  Theorems  5.1  and  6.1  establish  the 
asymptotic  and  spectral  exactness  of  the  estimator  S  in  (4.8). 


7.  A  posteriori  error  estimation  In  the  L_-norm 

In  the  present  section,  we  derive  a-posteriorl  estimators  for  |e.|  0 

IT(Q) 

To  this  end  we  consider  the  bilinear  form 


(7.1) 


Bj(u, v) 


Av  dx 


•  • 

on  Hj  x  Hg  where 


Hj  ■  Tu  «  L2(Q)  u(x1,x2,x3)  dXg  ■  0  a.e.  (Xj.Xg)  c  «|  , 


furnished  with  the  norm 
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|u|j  -  [J  lul2  dx) 


|U| 


L2(Q) 


and  where 


4  -  {v  .  H 


(Q) 


dv 


|Av|  <  •  ,  ■  0  on  R 

L2(Q)  an 


*} 


with  H  defined  5n  (i,  >) ,  furnished  with  the  nor* 
9 


|v|. 


( 


|Av|2  dxj 


1/2 


|Av! 


L‘(Q) 


We  note  that  locally  v  c  H2 (Q)  and  hence  is  well  defined.  Furthermore, 


it 


is  also  readily  seen  that  |«l2  is  a  norm  on 


Theorem  7.1.  The  bilinear  form  BjCu.v)  in  (7.1)  satisfies  (3.1)  and  (3.2) 
with  C  *  i  ■  1. 


Proof.  It  is  easy  to  see  that  (3.1)  holds  with  C  *  1.  We  will  now  estimate 


r  in  (3.2). 

• 

For  given  u  €  Hj,  define 

s  to 

be  the  solution  of 

(7.2) 

As  *  u 

in 

D. 

(7.3) 

s  *  0 

on 

r. 

(7.4) 

—  B  0 

an 

on 

R±  • 

2 

Since  u  e  L  (Q),  s  obviously  exists  and  is  uniquely  determined.  Define 


(7.5) 


Therefore 


.d/2 

:«s-Ij  . 

-d/2 


s(xrx2*x3)  3 * 


Bj(u,v)  «  I  uAz  dx  ■  |u|2,  |z|2  -  lu|j 
”*0 


and  (3.2)  follows  with  y  ■  1. 


□ 
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Me  turn  now  to  the  derivation  of  the  a-posteriori  estiaator  for  |e(q)|^ 

*  |e (q) |  „  Me  start  froa  the  characterization 
Lr(Q) 

(7.6)  B^e^qJ.v)  «  R^v) 


with  R^(*)  as  in  (4.16).  Using  Theorea  7.1  gives 

lRi  (Vj ) | 

ie.i  2  *  sup  rm - 

1  irfn}  v  2 

l  mj  V  L^(Q) 

■ 

where  the  supreaua  Is  taken  over  all  0  *  v  €  Hg  n  which  satisfy  (2. 13) 

* 

with  b  *  i.  To  estiaate  the  supreaua,  we  observe  that  any  v  e  Hg  can  be 
written  in  the  fora 


(7.7) 


'i'VVV  "  Z  ‘tk.V'VWfr] 


k.akl 


where  and  are  as  in  (6. 5)-(6. 8) .  Further,  we  find  that 


(7.8) 


where 


-4vi  *  Z  bik»  Vxfx2,*i.fr) 


k.akl 


bik.  ‘  aik.  Aik.  1  -  l’2’  «  "■ 


with  as  in  (6.8).  We  insert  (7.7)  in  Rj(Vj)  in  (4.12)  and  find  for 

1*2  with  q  *  2a  that 

W  “  r2(xl-X2){V2(xl*X2*d/2)  +  v2(xrX2’"d/2)  + 

JQ 

Z  42,2J  [  v2(xl'Vx3,L2jfr)dx3}‘,xldx2 

J*0  J-d/2 


where 


*2,2j  ’  -  a 
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Next,  using  (7.7)  yields 


W  -  2  I  “2M  "a  -  W 

k  .1 


where 


(7.9) 


T20£  “  0 


Vt 


*  pi 

T2ql  2  X  (4J+1)  4j( 

J=0  J-1 


z)  cos  (hrz)  dz 


where  q  *  2m  >  0.  Therefore 

W2*4  Z'oa)  T’i 

k  k 

where 


and  we  estimate 


°ti  ”  ^b2k«<A2k«*  <('1)  "  T2q^! 


4  S  Z  <b2k«)2  I  (1  -  *2,/  (W2  • 


Since  from  (7.7) 


(7.10) 

we  find 


,av>£(0,  *  s  I  « w 

k.m 


"e2'22 


sup 


W1' 


L&(Q)  v  |Av  |% 

Z  LZ(Q) 


5  d3,r2,22(u)  I  «-»>*  -  T2q<>2  (^’'4 


With  an  analogous  reasoning  for  i  ■  1,  we  have  shown 
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Theorem  7.2.  For  1  -  1,2 


|e.|  ?  S  VET  d3/2  |r  | 

1  L2(Q)  iq  1  L2(«) 


where 


Ei,  -  s  X  ((-1)t  -  Ti,t>2  ‘"it’ 


-4 


and  x2qi  ls  *n  t^.9),  and  for  1*1,  q  *  2a+l, 


■  pi 

Tlq«  '  \  Z  <4J*3)  i  4>»(z)  Sln  '"I*2'  d2' 

J*0  -1 


Remark  7. 1. 

1  -4 

For  q  *  0,  =  0  and  hence  EgQ  *  g  *  <(4)  *  1/180.  For  q  >  0, 

E,  ^  must  be  computed  numerically  in  general. 

8.  Error  estimation  for  laminated  materials. 

The  analysis  In  the  preceding  sections  carries  over  to  the  general 
problem  (1.1)  with  minor  modifications  which  we  will  describe.  To  underline 
the  analogy,  we  assume  that  a  basis  has  been  selected  so  that  its  span 

coincides  with  that  obtained  from  (2.6)-(2.9),  and  further,  that  for 

J  *  0,1,2,... 

(8.1)  *2J(Z)  *  *2J(_z)’  *2J+1(Z)  *  "  *2J+1(_Z)  * 
with 

(8.2)  b(z)  *j(z)*k(z)  dz  -  , 

J-l 

the  normalization  satisfied  by  L j ( z ) .  Then  we  obtain,  because  of  the  way 
span  <^>j>  was  defined,  that  for  1  *  1,2 
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(8.3) 


HjCv)  *  IVxi*x2){v*xi*x2*<i/2*  ±  vtej.Xg, -d/2)|  dXjdXg 

Jw 

ZV)  f  rl<xl,X2)  f  b^)*J  ‘WV  dx3dxldx2- 

<«n  J  it  J  _^a  /o 


J*0  « 


where 


corresponds  to  1  ■  1,  to  i 


v  e  S(q),  we  find  fro*  (8.2)  that 


»U  *  '  3  V1’  (2J“> 


2.  Since  R^v)  *  0  for  all 


OS  J5q 

1*1  if  J  is  odd, 
1*2  if  J  is  even. 


The  residuals  r^  are  defined  for  1  *  1,2  by 


(8.4) 


flu^f.q) 


ri*xl,x2*  *VX1,X2*  ”  dn  (xl*X2’d/2) 


Assuming  9  *  {«>,  we  reason  as  in  (4. 18)-(4. 20)  and  find 


(8.5) 


where 


rQ  S  f2r2  sup  (*i [ v)  )2  tiXjdXg 

V€H 

a.% 


♦jlv]  (Xj.Xg)  :> 


vlx^.d/Z)  ±  vlx^.-d/Z)  '  Y.  *1J  bm*jF?)v  dX3 

J*0  ^ -d/2 

frd/2  ,  1 1/2 

{  ‘|rf  |c)  (WV 


where  the  supreaua  is  taken  over 


M  :*  L 


-d/2  , 

(«;H1(-d/2,d/2)  n  |v  b^g-jv  dx3  *  0  a'*'  ^xl,x2^  €  W}  ’ 

J  /O 


Once  more  the  variational  problems  sup  ♦. [v]  adslt  unique  saxlalzers  v 

M  }  1 
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(8.9)  niqCXj.JCg)  *  ^  Clq  fCXj.Xg)  r^Xj.Xg)  . 

It  can  also  be  shown  that  exact  analogs  of  Theorems  5. 1  and  6. 1  hold  with 

suitably  modified  constants  D  .  We  shall,  however,  not  elaborate  since  the 

9 

details  are  completely  analogous. 

Romarh;  ft.  1 .  in  the  practically  important  case  that  a(z),  b(z)  are  piecewise 

constant  functions,  C,  in  (8.7)  can  be  computed  numerically  by  maximizing 

iq 

*[ v]  over  piecewise  polynomials  in  one  dimension  (of  sufficiently  high 
degree). 
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9.  A  staple  exaaole  of  a-posterlori  error  estimation 
We  consider 

(9.1)  Q«  (-1,1)  x  (-d/2,  d/2) 

and  select  f+  «  f  In  (1.1)  so  that 

u  -  u(x, y)  «  2  cos  |g  xj  cosh  fe  yj  . 

Then,  for  unifora  aodel  order  q  *  2a 

u(q)  -  £  UjUlL^  [f]  -  £  Xj  cos  [|  x]  Uy  [|^j 


J«0 


Then  the  vector  x  «  (xQ, 


j-0 


,  x  )  is  deterained  from 


(i1  *  4  1 1)  *  *  «• 


where 


and 


a  »  2*  sinh 


M  e 
(4  J  *  - 


(1,...,  1)  , 


[  4i4j  ®ij  -  [  4i  4j 

j_i 


dz 


Selecting  the  weight  function  f  »  1,  we  find 


and  the  estimator 


|e(q)|j  ^  «  a^2  cosh  -  xTej 


■  2 

^-inlnT  [“  -  a  Z  xj^jc1>)  '• 


1-1 


Using  a  coaputer  algebra  systea,  we  obtain 

2 

B2  f2 /|«(q)I j  f  -  1  +  d2  5-  ♦  0(d4) 
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where  a  is  listed  in  Table  9.1. 
<1 


B 

0 

2 

6 

8 

10 

12 

B 

240 

360 

936 

1768 

2856 

4200 

5800 

Table  9.1.  a^  is  the  asyaptotic  expansion  of  the  effectivity  index. 


Not  only  is  6  asyaptotically  exact  as  predicted  in  Theorea  5.1,  but  we 
observe  that  with  Q  =  0  and  5  =  a/2  we  have  k21  *  1  in  (5.3)  and 
in  (5.4) 

p  j2  2 

<*12>  *  1  ♦ - H - 

8{(2q+3 r  -  4} 

and  a  comparison  with  Table  9. 1  shows  that  for  q  i  2  this  bound  for 
k22  is  the  best  possible  one. 

Further,  it  is  verified  directly  that  Pj  -  cos  jg  xj  in  this  case, 
hence  #  in  (6.16)  is  equal  to  zero,  and  by  Theorea  6. 1  we  have  in  (5.3)  that 

(Aj)2  S  |  Dq  ,  1  «  1,2, 

i.e.  S  is  spectrally  exact  for  weight  functions  p  satisfying 

Q  -  0(qpd_p)  ,  0  5  p  <  1. 


Finally,  in  Table  9.2  we  present  the  asyaptotic  expansion  as  d  — » 0  of 


d3|r|2  /  |e(q)| 

LT(.u) 


Q)  ' 
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q 

d3  |r|2  /|e(q)|2 

l  («)  ir 

(Q) 

0 

180 

+ 

15(*d)2  / 

7  + 

0(d4) 

2 

630 

+ 

315(«d)2  / 

44  + 

0(d4) 

4 

2574 

+ 

1287 (sd) 2  / 

140  + 

0(d4) 

6 

6630 

9945(sd)2  / 

836  + 

0(d4) 

2 

Table  9.2.  Asymptotics  of  the  L  -residual  versus 
2 

the  L  -error  for  small  q. 


In  each  case  the  leading  tera  agrees  with  the  numerical  value  for 

obtained  from  Theorem  7.1  which  shows  the  asyaptotic  exactness  of  the 

estiaator  there  for  our  model  problem. 

Further,  in  the  unweighted  case  (l.e.  pal)  we  find  that  f  e  T_  with 

d2  f 

^  *  VX  \  where  A  is  the  first  eigenvalue  of - ~  in  (-1,1)  with 

dx 

boundary  conditions  u(±l )  «  0,  so  that  we  have  here 

f  «22  11/2 

*11  -  1  s  0  5  l1  +  r1  DqJ  *  *12 

with  as  in  (5.5),  i.e.  for  this  problem  the  estiaator  (4.8)  with  (4.25) 

is  asymptotically  and  spectrally  exact. 
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•  To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  BabuSka, Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742-2431. 


